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ABSTRACT 


The Hess panel clode was investigated as a procedure to predict the 
aerodynamic loading associated with propeller slipstream interference on 
the airframe. The slipstream was modeled as a variable onset flow to the 
lifting and nonlifting bodies treated by the code. Four sets of experi- 
mental data were used for comparisons with the code. The results indicate 
that the Hess code, in its present form, will give valid solutions for 
nonuniform onset flows which vary indirection only. The code, presently, 
gives incorrect solutions for flows with variations in velocity. Modifica 
tions to the code to correct this are discussed. 
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Influence matrix used In panel code solution 
Lift coefficient 

Normalized lift coefficient C. /C, 

^ ^MAX 

Propeller thrust coefficient. Reference (6) 

Propeller thrust coefficient, Reference (3) 

Unit normal vector on body surface 

Difference in static pressure between point on surface and 
free stream 

Slipstream dynamic pressure 
Free stream dynamic pressure 
Slipstream velocity 
Free stream velocity 
Angle of attack 

Panel source density used in panel code solution 
Slipstream dynamic pressure function used in Reference (5) 

Width of slipstream 


INTRODUCTION 


The renewed interest in propeller propulsion systems, stimulated by 
possible application to future transports, as well as the recognition of 
aerodynamic problems with existing twin engine aircraft, has created the 
demand for* analytical procedures to treat the combined slipstream/airframe 
flow field. A number of methods have been developed over the years to 
determine the flow solution for slipstream/wing combinations. Many of 
these were created within the last twenty years to support the interest 
' in deflected slipstream short-takeoff and landing vehicles. In all of 
these methods, however, the wing was treated as a lifting surface, i.e., 
having no geometric thickness. The effects of thick nonlifting bodies, 
such as the fuselage and nacelle, were not accounted for. It is well 
known that the system propulsive efficiency depends on the resultant flow 
about these bodiesj as well as the operation of the propeller itself. 

The most practical approach to analyzing the aerodynamic behavior of 
combined lifting and nonlifting bodies is to use a "panel code". This is 
a computer program which can calculate, with certain limitations, the flow 
field and corresponding aerodynamic loads about arbitrary bodies. A study 
was undertaken to investigate the prediction accuracy of one of these panel 
codes, the Hess code,^ compared to existing slipstream/wing aerodynamic test 
data. The operation of the code and the results of the investigation are 
reported in the following. 


THEORY 


Background 

The Hess code is a "panel method" computer code which calculates the 
linear potential flow about arbitrary three-dimensional lifting bodies. 

The geometry of the various components, which comprise‘the body, are input 
as corner coordinates of a set of trapezoidal panels which approximates the 
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surface of the component. Routines** are available which calculate these 
coordinates, given the basic component geometry and a paneling distribution 
scheme. 

Mathematically, the code solves the integral problem 

. n = -CT + /g3a/3n ds (1) 

When the surface is divided up into panels, on which the respective source 
densities are assumed to vary in a prescribed manner, the surface integration 
is performed piecewise on each panel, and equation (1) reduces to 





V 


( 2 ) 


This is an Nth order system of linear equations, for which there are numerous 
solution algorithms. The A., influence matrix is formed by determining the 
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normal component of velocity at control point i due to panel j with a unit 
source density. When there are lifting components, the influence of circula- 
tion equivalent unit dipole density is included from the respective lifting 
panels. After equation (2) is solved, the circulation equivalent dipole 
strengths are determined through application of the Kutta condition. The 
term V is called the onset flow and represents the free stream flow seen by 

po 
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the body. Normally, this term has the magnitude |V^| = 1. The vector 
components depend on whether angles of attack and sideslip are specified 


as part of the input. Equation (2) implies the is a constant, having ■ 
the same value for every panel on the body. This is not the case, however. 

The term V can take on different values for different panels, and accordingly, 

00 

can be used to represent a shear flow acting on the body. This applies to 
the case where it is required to represent the effects of a propeller slip- 
stream on the part of the body immersed in the slipstream. Equation (2) 
can be written in a more general form 


N 

l A.. 
j=l 




( 3 ) 


where V . indicates that the onset flow for panel (i) is used. An illustra- 
tion of this procedure is given in Figure 1. Figure 1(a) shows a paneled 
representation of a wing and nacelle with the conventional uniform onset 
flow. In Figure 1(b), the onset flow is varied for different streamwise 
panel strips to represent a propeller slipstream. Each of the respective 
onset flow strips in Figure 1(b) represents a V . vector in equation (3). 

The magnitude and angle of attack of the onset flow strip are determined 
from the axial and swirl flow components of the particular slipstream flow. 


Application of the Theory 

The purpose of the project was to test the Hess code, in its present 
form, as a procedure fpr calculating propeller slipstream induced airloads 
on fuselage/wing/nacelle configurations, A literature survey was performed 
to identify applicable experimental data for comparison. Four different 
documented experiments were selected, and the pertinent geometric and test 


condition data were coded for input. The results of the comparison cal- 
culations are presented later. 

As part of the project, a copy of the Hess code was provided by NASA 
for the necessary calculations. The code required some modification to 
make it compatible with the IBM based operating system of the university 
computer. During this phase, it became evident that gains could be made 
in reductions of the core storage requirements and execution times. These 
modifications are described in the following. 



CODE MODIFICATIONS 


The original version of the Hess code was developed on an IBM 370 
computer using the IBM FORTRAN IV language. The version which was made 
available by NASA for this project, had been modified to include the 
DERIV input geometry package coded in CDC FORTRAN. DERIV is described 
in Reference 2. This version was subsequently modified to run on IBM 
FORTRAN H, Optimization Level 2. The modifications consists of two 
parts; conversion of the CDC FORTRAN logic to equivalent IBM FORTRAN 
logic, and compression of core storage requirements and reduction in 
execution time. An additional modification was subsequently made to 
facilitate the input of the nonuniform onset flow velocity. This is 
discussed later in the report. The final version required only 268 kilo- 
bytes of memory in an overlay structure, and reduced execution time by a 
factor of two for a full one-thousand panel case. 

The remainder of this section is devoted to a detailed description 
of the logic and compression modifications made. It is intended for 
those who wish to make similar modifications. It assumes that the reader 
is familiar with the program structure and has a pertinent FORTRAN source 
listing. 

CDC Logic Conversion 

The modifications to the CDC logic involve four elements. The first 
of these is the use of a positive-indefinite variable for input array test 
purposes. The CDC octal code for this value is 17777000000000000000B. The 
equivalent IBM hexadecimal value is ZF7000000. In terms of real numbers, 
this code denotes a combination of a zero mantissa with an infinite exponent. 
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This value cannot be generated arithmetically or through Input conversloni 
therefore It serves Ideally as a test variable. 

Particular Input arrays are loaded with this value through the entry 
INDEF 1n subroutine ZERO prior to Input. After Input* the presence of this 
value is used to compress the Input data In subroutine SQUEEZE prior to 
output on a mass storage device. It is also used for test purposes In 
PANELIN. Here, the test variable nwst be integer. The test Is Invalid 
between two real variables because the ZF7000000 code Is converted to a 
floating point zero prior to performing the test. Also, In regard to this 
test In SQUEEZE and PANELIN, the masking part of the test used In the CDC 
code, can be removed with no difficulty. A questionable statement occurs 
In the DO 45 loop in PANELIN in this regard. A conversion Is made to a 
variable KHECK as part of an IF test. This statement does not work, as 
required, in IBM FORTRAN and must be removed. 

The second element concerns the use of the LOCF function throughout 
the code. Its purpose Is to provide the memory address of the argument. 

Its only essential use Is in the subroutine ZERO. Calls to ZERO and INDEF 
involve a string of memory from a conrion block whose length is detemnlned 
by specifying an Initial argument and a final argument In conjunction with 
LOCF. Subroutine ZERO and its associated calls were recoded so that the 
initial argument and the number of arguments of the common block string are 
used. A listing of the original and nwdified ZERO and CALL’S are given in 
Figure 2. 

The third element is the conversion from CDC specialized mass storage 
I/O (Input/output) to IBM random access I/O. This involves the subroutine 
READRS with entry WRITRS. The CDC system includes record keeping which must 
be explicitly coded for conversion. The OPENMS statement in NUED is replaced 
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by DEFINE flit. Also, the CDC I/O functions READMS and WRITMS are replaced 
by explicit coding in RANDIO. Because of the range of array sizes involved, 
the modified READERS subroutine stores the smaller arrays in core and writes 
the rest to mass storage using random access I/O. The storage and record 
keeping arrays are dimensioned for a maximum of five bodies (nonlifting) 
and five panels (lifting). The subroutine, its associated calls, and var- 
iable initializations are listed in Figure 3. 

The fourth element is a coding change in subroutines VFMNLF and VFMLFT 
necessary to run optimization level 2. This level is sensitive to the 
placement of ASSIGN and computed GO TO statements. Identical sequences 
of these are used in both subroutines. The coding was changed to use 
arithmetic IF statements, The old and modified coding are given in Figure 4. 

Core Storage Compression 
and Execution Time Reduction 

Figure 5 shows the overlay structure of the program with original 
significant common and/or dimension memory requirements in words. The 
important block is the 20,000 area in the subroutine COLSOL, which is the 
large A. ■ normal velocity matrix solver. This 20,000 block is used for 
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the core-resident part of the matrix which is being worked on by the 
solution algorithm. A large amount of mass storage I/O is required to 
shuttle pieces of the matrix to and from this block. The compression 
modification more than doubles the core-resident space for tho matrix. 

This in turn reduces the number of required disk I/O operations. For a 
one thousand panel solution on the university Amdahl V6 computer, the 
execution time was reduced by a factor of two. The modification is one 
of replacing all large memory blocks by a single labeled common block 
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dimensioned 42,900. The first 16,500 locations are used for labeled 
conmon blocks VRONSF, CTABLE, NORMAL, ATABLE, XYZAFT, and STABUv. To 
make this area available for COLSOL, the first 15,400 words are written 
on I/O unit 3 prior to the onset flows. The SIGMA array is not required 
until after COLSOL, The use of unit 3 for temporary storage is based on 
the desire to reduce the core required I/O buffers. In this regard, units 
12 and 15 have been removed from the program. Unit 12 was used for the 
second half of the source velocity matrix. Most disk packs can hold the 
entire matrix and there is no gain by using two I/O units. The entire 
matrix ^ss written on unit 11. Unit 15 was used in COLSOL to store the 
lower triangular part of the natrix solution for future use. In its 
present form, the code has no matrix solution for future use. The code, 
also has no provision for solving different right-hand-sides on separate 
runs, and hence, there is no requirement for maintaining this data on 
unit 15. Returning now to the temporary storage of the 15,400 words on 
unit 3, this is done in VFORM prior to writing the onset flows. For 
subsequent onset flow reads from unit 3, the data set is positioned by 
reading the known 15,400 word record back to its location in the common 
block. After COLSOL has been completed, the 15,400 word block is read 
back in VELPRS. The data set is then positioned for subsequent onset 
flow reads. With the use of the large common block, all array equivalencing 
is eliminated to take advantage of the optimizer. Care must be taken in 
PKUTTA where the logic requires PP and DELB to be equivalenced, or one 
substituted for the other. This is also the case in VELPRS for B and DD, 

The modified coding relating to the use of the large common block is given 
in Figure 6. 


11 


TEST CASES 


Four test cases were run. The requirements for each case were 
published data of experimentally determined spanwise lift distributions 
of a wing immersed in a slipstream flow. In two of the cases, the lift 
distributions were measured using spanwise wing segments with individual 
force transducers. In the other two cases, surface static pressure 
measurements were used. The published description of the test conditions 
varied from inadequate to fully adequate. The test cases and resulting 
comparisons with the Hess code are discussed in the following. 

Panel Modeling and Program Input 

The DERIV geometry routines were used for input and panel modeling. 

A general panel modeling scheme was applied to all of the test cases. 

The scheme was based on the 1100 maximum panel limit by the version of 
the code used. 

The propeller afterbodies, when present, were modeled as a duodecogonal 
cylinder. The cylinders were faired into the wing using straight-line taper 
in the vertical (Z) direction only. The afterbody chordwise panel density 
was made similar to that of the wing where the geometries overlapped. The 
entire spanwise section of the model occupied by the afterbody was assumed 
as nonlifting. 

The wing components were paneled according to the spanwise lift dis- 
tribution measurement scheme of the respective test models. Where in- 
dividual force measurement wing segements were used, the spanwise panel 
boundaries coincided with the segment boundaries. Where chordwise pressure 
distribution measurements were used, the spanwise panel boundaries were 
selected so that the panel control points coincided with the spanwise 


location of the pressure port rows. The chordwise panels were located 
according to the cos/sinh distribution function available in the DERIV 
procedure. Where possible, twenty chordwise panels on each of the upper 
and lower wing surfaces were used. Previous experience with the code 
indicated that this number was sufficient to give chordwise pressure 
distribution results of acceptable detail and accuracy. 

A modification to the code was made to improve the input of the 
nonuniform onset flow. In the original version, a delta velocity vector 
is input for each panel, once for the entire test case run. The final 
onset flow vector array used, is the vector sum of the uniform onset 
vector array (magnitude=l) and the input delta velocity vetor array. 

This means that when a particular nonuniform onset flow is desired, a 
delta input flow must be calculated which will produce this desired 
flow within the program. A second problem is that only one onset flow 
could be entered for each run. The modification solved these problems 
by shifting the delta velocity input to the only point in the program 
where the values are required. This point is in subroutine VFORM. Here, 
the final onset flows for the matrix solution are being formed for each 
of the input angles of attack. The input of delta velocity at this point 
allows a different nonuniform onset flow for each angle of attack, i.e. 
multiple nonuniform flow inputs per run. Also, the coding is altered so 
that the actual desired nonuniform flow is input and used as such, instead 
of the previous inte'^mediate delta flow. To avoid the inputting of a set 
of values for each of the 1000+ panels, the first and last index numbers 
for all of the panels in a particular streamwise strip are input with the 
respective nonuniform velocity vector values for that strip. This is the 
model shown in figure 1(b). 
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In three of the four test cases, usable slipstream data were avail- 
able. In the fourth case, the slipstream was calculated from given 
propeller geometric and operating data. The slipstream was input to 
the program as a set of nonuniform flow vectors as shown in Figure 1(b). 

The vectors were parallel to the chordwise plane, having an angle of 
attack only. For the calculated slipstream case, contraction effects 
were Ignored because of the relatively low loading of the propellers. 

For the available slipstream data cases, the data were taken well down- 
stream of the propeller plane where contraction is of no consequence. 
Accordingly, the slipstream input velocity vectors had no sideslip angle 
component . 

Test Case I , Reference (3) 

This test case was taken from Reference (3). The test model, shown 
in Figure 7, is a reflection plane configuration consisting of fuselage, 
wing and nacelle. The slipstream flow is generated by propellers. The 
wing was divided into eight segments, each providing direct force measure- 
ments of lift, drag and pitching moment. Two propellers having different 
pitch and chord distributions were used. The slipstream characteristics 
for the test conditions were inadequately described. There were sufficient 
propeller blade geometry and operating data to calculate the theoretical 
slipstream characteristics. 

The propeller performance method of Reference (4) was used for this 
purpose. As part of the solution procedure, the method determines the 
radial gradients of thrust and torque. These values were calculated for 
the conditions described in the reference, and then used to compute the 
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axial and swirl velocity components in the far wake by classic momentum 
analysis. 

The analytical model was paneled using seventeen (17) uniform span- 
wise strips. Each strip contained twenty (20) chordwise panels on the 
upper and lower surfaces respectively, The fuselage body was ignored, 
and replaced by a wing panel. The nacelle was paneled as a twelve (12) 
sided polygon with the longitudinal paneling concentrated near the lead- 
ing edge of the wing. The nacelle paneling increased in length toward 
the propeller and toward the trailing edge from this point. The complete 
spanwise section of the wing occupied by the nacelle was treated as non- 
lifting. The model paneling is shown in Figure 8, 

The test conditions calculated were for wing angles of attack at 0® 
and 10®, The propeller operating condition was specified in the reference 
as a thrust coefficient defined in terms of disk loading and slipstream 
dynamic pressure. A mid range value was selected for the slipstream character- 
istics. The comparative data were taken from Figures 24 and 42 in Reference 
(3). The comparisons are shown in Figure 9. The trends evident in the 
figure hold for the remaining test cases, The code generally over- 
predicted the lift magnitudes while giving lift distributions similar to 
the experimental. The errors in the lift magnitudes increased with test 
condition angle of attack. This is not unexpected as the code does not 
account for loss of lift due to separation and stall, The discrepancies 
in lift distribution are due in part to the use of calculated slipstream 
characteristics rather than actual measured data. 
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Test Case It» Reference (5) 


The test model is shown in Figure 10, The confisuration is of an 
infinite aspect ratio wing, The slipstream is produced by a propeller 
with drive system located upstream. The lift distribution is detemined 
by direct force measurements from the individual spanwise segments* The 
characteristics of the slipstream are well documented In tenns of both 
axial and swirl components. 

The iimthematical model was paneled with twenty (20) spanwise segments. 

Each segment contained twenty (20) chordwise panels on the upper and lower 
surfaces respectively. The spanwise paneling distribution was the same as 
the experimental model except that the narrow width strips were continued 
across the slipstream through the spanwise wing center, The ntodel panel 
layout Is given in Figure 11, 

Six different test conditions were calculated. All were at an angle 
of attack of and included a no^slipstream condition. The spanwise 
lift coefficient distribution used in the reference was based on a function 
of both the free stream and the slipstream dynamic pressures. No correct- 
ing factor could be developed from the reference which would adjust the 
calculated results for this. Therefore, a normalised lift coefficient was 
used for comparison. All experimental data were normalUed by the maximum value 
of the experimental lift coefficient for the no-slipstream flow case. The 
cal^railated data were correspondingly normalised by the calculated no-slip- 
stream maximum lift coefficient. Figure 12(a) shows the no-slipstream case 
in terms of the respective nontialised lift coefficients. Figures 12(b) and 
12(c) give the comparisons for two of the five remaining test conditions. 

The 0 values listed are the aforementioned dynamic pressure functions* 
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Increasing c? indicate increased slipstream velocities. The computer 
code again is approximately modeling the distribution but overpredicts 
the magnitude within the slipstream. The overprediction increased with 
increasing slipstream velocity, 

Test Case III, Reference (6) 

The experimental model, shown in Figure 13, is similar to the model 
of Test Case I. It is a reflection plane model consisting of a wing, 
nacelle and propeller. The propeller span extends beyond the wing tip. 

The lift distribution was measured by surface static pressure distribution. 
The slipstream characteristics were adequately described in terms of both 
the axial and swirl components. The wing section was symmetrical and the 
test conditions included negative as well as positive angles of attack. 
Different propeller thrust levels as well as no-slipstream were run. 

The mathematical model was paneled so that the spanwise segment centers 
aligned with the chordwise rows of static pressure ports. Because of the 
number of panels used, the chordwise paneling was reduced to seventeen (17) 
instead of twenty (20). The nacelle was paneled similar to Test Case I. 

The model paneling is shown in Figure 14. 

The comparison results are given in Figure 15. Figure 15(a) is the 
no-slipstream condition. Here* the code underpredicts the lift magnitudes. 
This underprediction also occurred with test case IV, and is discussed 
there. In Figure 15(b), at a low thrust coefficient, there is good agree- 
ment. However, in Figure 15(c) at a higher thrust coefficient, the mag- 
nitude overprediction occurs again. The distribution patterns are predicted 
quite well. 


Test Case IV, Reference (7) 


The experimental model is an infinite aspect ratio wing as in 
Test Case II. The model is a wing only. The slipstream was created 
by an axial flow blower located upstream of the wing, Figure 16. The 
slipstream contained only an axial component. There is no angle of 
attack variation due to swirl. The lift distribution was determined 
from surface static pressure ports. The test conditions consisted of 
variations in the axial flow velocity. 

The analytical model was paneled so that the spanwise segment centers 
aligned with the chordwise pressure port rows. Twenty (20) spanwise strips 
and twenty (20) chordwise panels were used. The model is shown in Figure 
17. Only the spanwise panel boundaries are shown. An illustration in* 
eluding the chordwise boundaries was not of sufficient clarity due to the 
high line density. 

The comparison results are shown in Figure 18. Again the code pre- 
dicts higher magnitudes in the slipstream region. Outside of the slip- 
stream, the code underpredicts the lift distribution. As with the propeller 
tests, the prediction error in the slipstream region increases with increasing 
slipstream velocity. 

Figures 19 and 20 give experimental and calculated chordwise pressure 
distributions respectively for spanwise locations inside and outside of the 
slipstream. The data are for the flow condition of Figure 18(b). It is 
repeated here that this test case is for a jet flow, which has no angle of 
attack variation, unlike the previous propeller test cases. Figure 19 shows 
the expected increase in pressure magnitude within the slipstream. It is 
noted that the pressure approaches free stream static at the trailing edge 
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regardless of the existence of the slipstream. The corresponding cal- 
culated results are given in Figure 20. There are two distinct variations 
from the experimental. First the slipstream distribution does not return 
to free stream at the trailing edge. A comparison of the two distributions 
in Figure 20 showed that the calculated airfoil chordwise velocity distribu- 
tion within the slipstream was an exact magnification of the velocity dis- 
tribution outside of the slipstream by a factor of the jet velocity ratio 
(V. ^/V = 2.5). The effect of the nonuniform onset flow input was that of 
a pure magnification. The calculated pressure distributions result from 
the set of velocities that are the solution of Hess code. Pressures are 
not dealt with specifically. If the inside slipstream distribution were 
corrected by subtracting the total pressure change due to the jet velocity, 
the trailing edge pressure would fall well below free stream static. In 
addition, there would still be a large difference at the leading edge 
between the calculated and experimental distributions. The conclusion 
offered is that the code does not calculate the correct lifting body chord- 
wise distribution for nonuniform onset flows. There is a circulation change 
occurring experimentally that is not accounted for as the code is presently 
formulated. 

The second variation between experiment and calculation is in the 
airfoil chordwise distributions outside of the slipstream. These are shown 
in Figure 21. There is considerably more disagreement for a pure airfoil 
problem than has been previously experienced by the author. Figure 21 shows 
the source of the underprediction seen in Figure 18. A possible explanation 
lies with the paneling of the airfoil used in the test. The airfoil is a 
NACA 0009, which has a smaller leading edge radius than the other airfoils 
used in the test cases. It is possible that the panel density should have 
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been increased in the leading edge region to obtain a better representa- 
tion of the rapidly varying noise curvature. For the underprediction in 
Test Case III, seventeen (17) panels were used for the airfoil instead of the 
selected twenty (20). Here again it is possible that the airfoil was 
under paneled for the desired accuracy. In both cases, the 1100 panel 
limit did not allow increases in the paneling, and the question of the cause 
of these variations is still unresolved. 
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SUMMARY 


Comparisons with the test cases show that the Hess code gives 
acceptable representations of the variation in lift induced by a 
propeller slipstream. The code tends to predict higher lift mag- 
nitudes than measured for increased angles of attack and slipstream 
velocities. These errors are due to the lack of a formulation within 
the code to treat the problem of a lifting body with jet flow. The 
code appears to give valid results for nonuniform onset flows where 
only the angle of attack varies, i.e. the velocity magnitude remains 
unity. 

If a modification to the Hess code can be made which will work the 
jet flow solution, then the code will be suitable for the slipstream/ 
airframe problem. A possible solution to this problem is to utilize 
singularity distributions to model the slipstream instead of the onset 
flow. The strength and distribution functions of the singularities are 
determined beforehand to produce the desired slipstream characteristics 
of the isolated propeller. These known singularities are then included 
into the Hess code formulations for the induced velocity influence co- 
efficients. The nature and geometric location of these singularities 
is only partly evident at this time. The axial flow component can be 
represented by a radial distribution of ring sources at the propeller 
disk. This is equivalent to semi-infinite co-axial vortex tubes which 
is the classic means of representing the axial flow of the propeller 
wake. The swirl component is considerably more difficult. The singular- 
ity surfaces, which generate this flow, will penetrate the body geometry 
and result in numerical difficulties when determining the induced flow 
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components. The modeling of this flow can be done either with vortices 
or doublets; however, there is insufficient experimental data to provide 
proper guidance in this matter. 

In summary, the results of this investigation support the use of Hess 
code to study propeller propulsion system integration problems. Care must 
be exercised at this stage when considering the validity of the code’s 
predictions for high lift and high thrust conditions. Improvements in the 
prediction accuracy seem possible by replacing the variable onset flow with 
singularity distributions. This, however, will require an experimental 
program to determine the most practical modeling of the slipstream swirl 
component. 
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(a) Uniform onset flow 



l-Vo, =1“^ 

(b) Variable onset flow 


Figure 1. Use of variable onset flow to represent 
propeller slipstream. 
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ORIGINAL PASS IS 
OF POOR QUALITY 


SUBROUTINE ZERO (A.B) 

DIMENSION A(l) 

STUFFsO.O 

10 N»LOCF(B)-LOCF(A)+1 

DO 20 1=1, N 

20 A( I )»STUFF.ANO. (LOCF(A( I ) ) .OR. 77770000000000000000B) 

RETURN 
ENTRY INDEF 

STUFF=17770000000000777777B 

60 TO 10 

END 

(a) original ZERO 


LOCATION 


CALL ZERO (DELOSX(l). DELOSZ(UOO) NUED 

CALL ZERO (LENGTH (1), LENGTH (199)) NUED 

CALL ZER0(XS(1),DZ(49)) BOOYIN 

CALL ZERO(BAREF,NYIG) BODYIN 

IF(BTHI.NE.l) CALL ZER0(VY(1) .VY(NVYP)) BODYIN 

CALL ZER0(XX(1),XZ(NXX,NYY)) B0D6E0M 

CALL ZERO(PAREF,ROTG) PANELIN 

CALL ZERO ( CHORD (1), CHORD (20)) PANELIN 

CALL INDEF(X0FF,Z0FF(200)) READIN 

CALL INDEF(VX(1),JYIG(50)) BODYIN 

CALL INDEF(PXL(1),VS(41)) PANELIN 

CALL INDEF A(LVX),A(LXX-1)) VORBOO 

CALL INDEF( ,A) VORPANS 


(b) original CALL'S 


Figure 2, Original and modified subroutine ZERO and CALL'S 
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10 

20 


(c) 


(d) 


Figure 2 


ORIGU'^HL U 

QH POOR QUAUTY 


SUBROUTINE ZERO(A.N) 

DIMENSION A(l) 

DATA PINDEF ,HEXZR0/Z7F000000 , ZOOOOOOOO/ 

STUFF=HEXZRO 

DO 20 1=1 ,N 

A(I)=STUFF 

RETURN 

ENTRY INDEF(A,N) 

STUFF*PINDEF 
GO TO 10 
END 

modified ZERO 


LOCATION 


CALL ZERO(DEL0SX(l),3300) 

NUED 

CALL ZER0(L£NGTH(1,1),30) 

NUED 

CALL ZERO (XS(1), 2008) 

BODY IN 

CALL ZER0(BAREF,29) 

BODY IN 

I F( BTH I . NE . 1 ) CALL ZERO (VY(l).NVYP) 

BODY IN 

CALL ZER0(XX{1),NXX+2*NXX*NYY) 

BODGEOM 

CALL ZER0(PAREF,66) 

PANEL IN 

CALL ZERO ( CHORD (1), 20) 

PANELIN 

CALL INDEF(X0FF,600) 

READIN 

CALL INDEF(VX(1),702) 

BODY IN 

CALL IN0EF(PXL(1),2202) 

PANELIN 

CALL INDEF(A(LVX),LXX-LVX) 

VORBOD 

CALL INDEF(A,1) 

VORPANS 

modified CALL'S 



. Concluded 


26 


XXX x: X X 


ORIGINAL PAGtt feS 

OF POOR QUALITY. 


COMMON/RANF/ INDEX(200) »LENGTH(199) , 

GENOA! , VORTEX ,C0NPTS .NORMAL , 
CONFUN .VORTIN .PLANFM .TWICAM . 
THICKS .CPOINT .ELEVEN .TWELVE 


CALL OPENMS( 19, INDEX, 200,0) 

(a) COMMON and initialization statements in NUED 

SUBROUTINE READRS(A,N,K,IC) 

COMMON/RANF/ INDEX ( 200) .LENGTH (199) , 

GENDAT .VORTEX .CONPTS .NORMAL , 

CON FUN . VORT I N . PLANFM .TWICAM. 

THICKS. CPOINT. ELEVEN .TWELVE 
EQUIVALENCE ( CON PNT. CONPTS) 

DIMENSION A(N) 

DATA LUN/19/ 

CALLMHREAD 

KEY=K 

IF(KEY.LE.12) KEY=12*-(IC-I)+K 
LEN=N 

IF(LEN.EQ.O) LEN 'LENGTH (KEY) 

IF(LEN.EQ.O) RETURN 

IF LENGTH ( KEY ) . EQ . 0 ) LENGTH ( KEY ) *LEN 

IF CALL.EQ.4HREAD CALL READMS(LUN.A.LEN.KEY) 

IF(CALL,EQ.4HWRIT) CALL WRITMS(LUN.A.LEN.KEY.IREWR) 

RETURN 

ENTRY WRITRS 

CALL=4HWRIT 

IREWR=0 

IF(N.EQ.O) IREWR=-1 

GO TO 10 

END 

(b) original READRS 


Figure 3. Original and modified READRS and CALL'S 


LOCATION 


POOR ^ 


CALL WR ITRS ( XOFF , SQUEEZE ( XO FF , 600 , 1 ) , ELEVEN , 1 ) 

CALL WRITRS(A,LEN,PLANFM,IBOO) 

CALL WRITRS(JXIG,SQUEEZE{JXir7.350,l),CP0INT,IB0D) 

CALL WR ITRS VX , SQUEEZE ( VX , 352 , 1 ) , VORTI N , I BOO ) 

CALL WRITRS(BAREF,29,GENDAT,IB00) 

CALL WRITRS(PXL,140,PLANFM,NB0D+IPAN) 

CALL WR ITRS ( PXC , SQUEEZE ( PXC , 1 020 , 1 ) ,TW I CAM ,NBOO+I PAN ) 
CALL WRITRS{X0CT,SQUEEZE(X0CT,960,1),THICKS,NB00+IPAN) 
CALL WRITRS (VC ,82 ,V0RTIN ,NBOO+IPAN) 

CALL WR ITRS ( PAREF , 66 , GENDAT ,NB0D+I PAN ) 

CALL READRS(BAREF,0, GENDAT, I) 

CALL READRS(A(LXX),0,PLANFM,I) 

CALL READRS(A,0,V0RTIN,I) 

CALL READRS(A(LXIG),0,CP0INT,I) 

CALL READRS ( PAREF, 0, GENDAT, I) 

CALL READRS (PXL,0,PLANFM, I) 

CALL READRS(VC,0,V0RiIN,I) 

CALL READRS (Z(KPX),0,TWICAM, I) 

CALL READRS(A(KXT),0,THICKS,I) 

CALL READRS ( A ( KXOFF ) . 0 , ELEVEN , 1 ) 

CALL READRS ( PAREF ,0. GENDAT, II) 


(c) original CALL'S 


COMMON /RANF/ I NDEX (6,5) , LENGTH (6,5), STORE (317 ,5 ) , 
X L0CINT(4),L0CFIN(4).IDACSS 


DEFINE FILE 19(384, 100, U,IDACSS) 
IDACSS*! 

L0CINT(1)=1 

L0CINT(2)=30 

L0CINT(3)=170 

L0CINT(4)=252 

L0CFIN(1)=29 

L0CFIN(2)=169 

L0CFIN(3)=251 

L0CFIN(4)=317 


(d) modified COMMON and initialization statements in NUED 


Figure 3, Continued 


READ IN 

BODY IN 

BODY IN 

BODY IN 

BODY IN 

PANELIN 

PANEL IN 

PANELIN 

PANELIN 

PANELIN 

VORBOD 

VORBOD 

VORBOD 

VORBOD 

VORPANS 

VORPANS 

VORPANS 

VORPANS 

VORPANS 

VORPANS 

BPCOEF 
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ORidll^AL PASS ^ 
OF POOR QUAUTY 


SUBROUTINE READRS(A,N,K JC) 

COMMON/RANF/ INDEX{6,5) ,LENGTH(6,5) ,ST0RE(317,5) . 
X L0CINT{4),L0CFIN(4),IDACSS 

DIMENSION A(N) 

DATA LUN/19/ .READ .WRIT/4HREAD .4HWRIT/ 

CALL=READ 

10 IF(K.GT.6) 60 TO 50 

LEN=N 

IF(LEN.EQ.O) LEN=LEN6TH{K,IC) 

IF(LEN.EQ.O) RETURN 

IF(LENGTH(K,IC).EQ.O) LENGTH (K,IC)=LEN 
IF(CALL.EQ.READ) 60 TO 25 
IF(IDACSS.NE.l) IDACSS=NXTWRT 
INDEX (K,IC)=IDACSS 
DO 15 I=1,LEN.100 
CALL WRITMS{LUN,A(I).IDACSS) 

15 CONTINUE 

NXTWRT=IDACSS 

RETURN 

25 IDACSS^ INDEX (K,IC) 

DO 30 I=1,LEN,100 

CALL READMS(LUN,A(I).IDACSS) 

30 CONTINUE 

RETURN 

50 II=L0ClNT(K-6) 

IF=L0CFIN(K-6) 

I F (CALL. EQ. READ) 60 TO 65 
DO 60 I=II,IF 
ST0RE(I,IC)=A(I-II+1) 

60 CONTINUE 

RETURN 

65 DO 70 I=II,IF 

A(I-II+1)=ST0RE(I.IC) 

70 CONTINUE 

RETURN 

ENTRY WRITRS(A.N.K.IC) 

CALL=WRIT 
60 TO 10 
END 


t (e) modified READRS 

I 

* 

I . ■ 


Figure 3. Continued 
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SUBROUTINE RANDIO 
DIMENSION A(IOO) 

ENTRY READMS(LUN,A,/IDACSS/) 

READ(LUN'IDACSS)A 

RETURN 

ENTRY WRITMS(LUN,A,/IDACSS/) 

WRITE(LUN'IDACSS)A 

RETURN 

END 


(f) replacement subroutine for CDC READMS and WRITMS functions 


LOCATION 


CALL WR ITRS ( XOFF . SQUEEZ { XOFF , 600 , 1 ) , 1 , 1 ) 

READIN 

CALL WRITRS(A,LEN,2,IB0D) 

BODY IN 

CALL WRITRS(JXIG,SQUEEZ(JXIG,350,1),3,IBOO) 

BOOYIN 

CALL WRITRS ( VX , SQUEEZ ( VX , 352 , 1 ) , 4 , 1 BOO ) 

BODY IN 

CALL WRITRS(BAREF,29,7,IB00) 

BODY IN 

CALL WRITRS(PXL,140,8,IPAN) 

PANEL IN 

CALL WRITRS(PXC,SQUEEZ(PXC, 1020,1) ,5, IRAN) 

PANEL IN 

CALL WR ITRS ( XOCT , SQUEE Z{ XOCT , 960 , 1 ) , 6 . 1 P AN ) 

PANEL IN 

CALL WRITRS(VC,82,9,IPAN) 

PANELIN 

CALL WRITRS(PAREF,66,10,IPAN) 

PANELIN 

CALL READRS(BAREF,0,7,I) 

VORBOD 

CALL READRS(A(LXX) ,0,2,1) 

VORBOD 

CALL READRS(A,0,4,I) 

VORBOD 

CALL READRS (A, LXIG) ,0,3,1) 

VORBOD 

CALL READRS(PAREF,0,10,II) 

VORPANS 

CALL READRS(PXL,0,8,II) 

VORPANS 

CALL READRS (VC ,0,9, II) 

VORPANS 

CALL READRS(A(KPX), 0,5,11) 

VORPANS 

CALL READRS (A(KXT), 0,6,11) 

VORPANS 

CALL READRS(A(KX0FF) ,0,1,1) 

VORPANS 

CALL READRS (PAREF, 0,10, I) 

BPCOEF 


(g) modified CALL'S 


Figure 3. Concluded 


ORIGINAL PAGE IS 

OF POOR QUAUTV 


IF ( NSYM - 1 ) 58,52,54 
52 ASSIGN 910 TO 119 
GO TO 60 

54 ASSIGN 920 TO 119 
GO TO 60 

58 ASSIGN 2000 TO 119 


/ 

/ 

DO 1700 12 = 1, LOOP / 

IF ( 12 .EQ.LOOP ) GO TO 119, ( 2000, 910, 920) 

GO TO (1000, 910, 920, 910 ), 12 

(a) original 


DO 1700 12 = 1, LOOP 
IF( 12. EQ.LOOP) GO TO 880 
IF( 12-2) 1000,910,870 
870 IF(I2-4)920,910,880 

880 IF(NSYM-1)2000.910,920 

(b) modified 


Figure 4. Original and nradified coding in subroutines VFMNLF 
and VFMLFT 


Of . POOR 



Figures. Overlay structure. Numbers denote size 
of DIMENSION and COMMON respectively 
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X X X X X X XXX XXX 


nu^in 


COMMON / VARONSF / OELOSX(UOO), DELOSY(llOO) , DELOSZ(llOO) 
COMMON / CTABLE / XC(UOO), YC(llOO), ZC(llOO) 

COMMON / NORMAL / XN(llOO), YN(llOO), ZN(llOO) 

COMMON / ATABLE / A(llOO) 

COMMON / XYZAGT/ XNGT(llOO), YNGT(llOO), ZNGT(llOO) ,AGT(1100) 
COMMON / STABLE / SICiMA(UOO) 

(a) original COMMON blocks 


COMMON / BLANk/DEL0SX(1100),DEL0SY(1100,KDEL0SZ(1100), 
XC(1100).YC(1100),ZC(110G), 
XN(1100),YN(1100),ZN(1100). 

A(llOO), 

XNGT(llOO) .YNGT(llOO) ,ZNGTtll0(J.AGT{1100) , 
SIGMA(llOO), 

DUMMY (26400) 


(b) modified COMMON block 


COMT^N /BLANK/DELOSX (1100), OELOSY (1100) .DELOSZ (1100), 
DUMMY (13200), 

XOFF(2O0.YOFF(2OO) ,Z0FF(200) , 

DUMMY2( 25800) 

(c) use in READIN and GPRINT 


COMMON /BLANK/DUMMYK 16500), 

XS(41) ,TYS(21 ,41) ,RZS(21 ,41) ,XMC(49) ,YM(49) ,ZM(49) , 
DY(49),DZ(49),VX(301).VY(51), 

JXIG(300) ,JYIG(50 ,DUMMY2(100) ,A(23590) 

(d) use in BODYIN and BOOPRT 


Figure 6. Modification to large array storage requirements 


33 


xxxxx xxxxxxxxx xxxx 






COMMON /BLANK/DUMMYK 16500), 

PXL(20) ,PYL(20) ,PZL(20) ,CH0R0(20) ,PXT(20) ,PYT(20) , 
PZT(20) ,PXC(30) ,PET(30) .CAMBER (30, 30) ,TWIST(30) , 
TWPCR(30) ,X0CT(30) ,ETAT(30) ,THICK(30.30) . 
VC(41),VS(41), 

DUMMY2(24198) 

(e) use in PANELIN and PANPRT 


COMMON /BLANK/DUMMY (16500) ,A( 26400) 

(f) use In VORBOD and VORPAN 

COMMON / BLANK / DELOSX(llOO) , DELOSY(llOO) , DELOSZ(llOO) , 



DUMMYl ( 13200 ), 




vxs 

(1100 ) , 

VYS 

(1100 ) 


VZS 

(1100 ) , 

DIPLEX 

1100 ) 


DIPLEY (1100 ) , 

DIPLEZ 

(1100 ) 


VXL 

(1100 ) , 

VYL 

(1100 ) 


VZL 

(1100 ) . 

VXR 

(1100 ) 


VYR 

(1100 ) , 

VZR 

(1100 ) 


VI 

(1100 ) . 

VJ 

(1100 ) 


VK 

(1100 ) , 

DUWiY2 

(9900 ) 


se in VFORM 






COMMON /BLANK/DUMMYl ( 3300 ) , 

XC(1100),YC(1100).ZC(1100), 

DUMMY2(9900), 

XIJ(llOO) .YlJdlOO) .ZIJ(llOO) , 
DUMMY3( 23100) 

(h) use in VFMNLF 


Figure 6. Continued 
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ORIGINAL PAGK B 
OF POOR QUALITY 


COIflON /BLANK/DUMMY 1 ( 3300 ) , 

xc(noo),YC(iioo),zc(iioo)c 

DUMMY2(9900), 

XU (1100 ) , YU (1100 ) 

ZU (1100 , VXL (1100 

VYL 1100 , VZL 1100 

VXR (UOO ) , VYR (1100 

VZR (1100 ) , VXTRAX 1100 ) 

VXTRAY (1100 ) , VXTRAZ (1100 ) 

DUMMY 3 (13200) 

(i) use In VFMLFT 


COMMON / BLANK / DELOSX(llOO) , DELOSY(llOO) , DELOSZ(llOO) , 
DUMMYl (13200), 

VSUMX (1100 ) , VSUMY (1100 ) , 

VSUMZ (1100 ) , VDIFX 1100 ) , 

VDIFY 1100 ) , VDIF2 (1100 ) , 

VLASTX 1100 ) , VLASTY 1100 ) , 

VLAST2 (1100 ) , VNEXTX (1100 ) , 

VNEXTY (1100 ) , VNEXTZ (1100 , 

ONSETX (1100 ) , ONSETY (1100 ) , 

ONSETZ (1100 ) , 

VXL 1100) , VYL (1100 ) , 

VZL ( 1100) . VXR 1100 ) , 

VYR ( 1100) , VZR (1100 ) , 

DUMMY2(3300) 

(j) use in PSWISE 


COMMON /BLANK/DELOSXdlOO) ,DEL0SY(1100) ,DEL0SZ(1100) , 
XC(1100),YC(1100),ZC(1100), 
XN(1100),YN(1100),ZN(1100), 

A(llOO), 

XNGT(llQO) ,YNGT(1100) ,ZNGT(1100) ,AGT(1100) , 
SIGMA(llOO), 


vx 

(1100 ) . 

VY (1100 ) , 

vz 

(1100 ) , 

RHSIDE (1100 ) , 

AU 

1100 ) , 

DIPLEX 

(1100) . 

DIPLEY (1100) , 

DIPLEZ 

(1100 , 

DUMMY (17600) 



(k) use in AFORM 
Figure 6. Continued 
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XX >c X X X X X>C>C>C ^ X X X xxxxxxxxxxxx 


OF POOR 


,li» 


COMMON /BLANK/A (42900) 

(1) use in COLSOL 


COMMON /BLANK/OUMMYl (15400), SI6MA(1100), 


FLOWX 

( 40 , 

FLOWY 

( 40 ) 

1 

FLOWZ 

{ 40 ) 

VELX 

( 40 ) 

1 

VELY 

( 40 ( 

VELZ 

( 40 ) 


SUMX 

( 40 ) 

SUMY 

( 40 ) 

• 

SUMZ 

( 40 ) 

VXINF 

( 40 ) 

t 

VYINF 

( 40 ) 

VZINF 

( 40 ) 


VXFLOW 

i 40, 20 ) , 

VYFLOW 

( 40, 20 ) , 

VZFLOW 

{ 40, 20 ) , 

DVSQ 

20 ) 

» 

TERMl 

( 20 ) 

S 

( 20, 20 , 

PS 

( 20 ) , 

P 

( 


DELB 

• 

o 

CM 

D 

( 20. 

20 ) . 

EMINF 

f 2 ) . 




EM 

( 2 . 20 ) . 

qiNF 

20. 20 ) . 

0(20,20,20) 
DUMMY2 (14178) 





(m) use in PKUTTA 


COfflON /BLANK/DUMMYl (6600) , 

XN(1100),YN(1100),ZN(1100), 

DUMMY2(1100), 

XNGT(llOO) ,YNGT(1100) ,ZNGT(1100) ,AGT(1100) , 
SIGMA(llOO), 

COMSIG(llOO), VX(llOO), VY(llOO), VZ(llOO), 
VEL(llOO), DCX(llOO), DCY(llOO), DCZ(llOO), 
VN(llOO), PCOEF(llOO), VSUM(llOO), 
ONSETX(llOO) ,ONSETY(1100) ,0NSETZ(1100) , 
caowx(lioo) ,CFL0WY(1100) ,CFLOWZ(llOO) , 
VXS(llOO) .VYSdlOO) ,VZS(1100) , 

UNIFMXdlOO) ,UNIFMY(1100) ,UNIFMZ(1100) , 
USIGMA(llOO) 

(n) use VELPRS 


Figure 6. Continued 
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COMMON /BLANK/OUMMYl ( 3300 ) , 

XC(1100),YC(1100),ZC(1100), 

XN(1100),YN(1100),ZN(U00), 

A(llOO), 

DUMMY2(1100), 

XNGT(llOO) JNGT(nOO) ,ZN6T(1100) ,AGT( 1100 ) , 
COMSIG(llOO), VX(llOO), VY(llOO), VZ(llOO), 
VEL(llOO), DCX(llOO), OCYdlOO), DCZ(llOO). 
VN(llOO), PCOEF(llOO), VSUM(llOO), 

MXI(200), MYJ(200), MZK(200), 

FSTRPX(200), FSTRPY(200), FSTRPZ(200), 
SECM0X(50), SECM0Y(50), SECM0Z(50). 
FSECX(50), FSECY(50), FSECZ(50), 

X^^■0T, YMTOT, ZMTOT, XFTOT, YFTOT, ZFTOT. 
DUMMY3( 12794) 


(o) use in PRINT 


X 

X 

X 

X 

X 

X 

X 

X 


COMMON /BLANK/DUMMY1(3300), 

XC(llOO), YC(llOO), ZC(llOO), 

DUMMY 2 (22000), 

MXI(200), MYJ(200), MZK(200), 

FSTRPX(200), FSTRPY(200), FSTRPZ(200), 
$ECM0X(50), SECM0Y(50), SECM0Z(50), 
FSECX(50). FSECY(50), FSECZ(50), 

XMTOT, YMTOT, ZMTOT, XFTOT. YFTOT, ZFTOT, 
DUMMY3( 12794) 


(p) use in BPCOEF 


Figure 6. Concluded 
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Figure 7. Reflection plane test 
model from Reference (3). 
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(a) Top view 


Figure 8. Panel layout for the test model of Reference (3). 
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(e) AiH‘o11 


'6 8, Concluded 







ORIGINAL PAGE 18 
OF POOR QUALITY, 


Experimental 



Figure 9. Comparison with Reference (3), 
calculated versus experimental. 

^a) Reference Figure 24, a=0°, Cy = 0*6 

's 


41 




OF POOR 





Figure 9. continued 

(b) Reference Figure 24, a=10“,CT =0.6 

I c 
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SPAN 

Figure 9. concluded 


(d) Reference Figure 42 , a=10°, C^. = 0.5 



ORIGINAL PAGE fS 

OF POOR QUALITY 



Figure 10. Infinite aspect ratio 
test mooel from Reference (5). 
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(a) Top view 



Figure 11. Panel layout for the test model of Reference (5) 



ORkGlMAl. PAGE S3 

OF POOR QUALllY 


— Experimental 
12 o Calculated 



-16 -12 -8 -4 0 4 8 12 16 

Span 


Figure 12. Comparison with Reference (5), 
calculated versus experimental. 

(a) Reference Figure 10a, a=6i“, a = 0.0 
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Figure 12. Continued 

(b) Reference Figure 10a, «=6i°, a = 0.26 
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Figure 12. Concluded 

(c) Reference Figure 10b, a=6i°, o = 0.49 


If A 


POOR OUi^LITV 


a 



Figure 13. Reflection plane test 
model from Reference (6). 




(a) Top view 


Figure 14. Panel layout for the 
test model of Reference (6). 
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(b) Side view 



« 


^ 5 ' POOR 


EKperlmental 



Figure 15. Comparison with Reference (6), 
calculated versus experimental. 

(a) Reference tables 5 and 6, a«5®,10®, C^«0 


\ 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



SPAN 

Figure 15. Continued 


(b) Reference tables 7-11, 


a=-10°, -5% 0 \ 5 \ 


10°; C^=0.36 
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omGlNftL Pft® S 
OF POOR QUM-ITY 



Figure 15. Concluded 


(c) Reference table 12-16, a=-10°, - 5 °, 



C^=0.64 
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Figure 16. Infinite aspect ratio 

test model from Reference (7). 
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PAGE JS 



(a) Top view 


Figure 17. Spanwise panel layout 
for the test model of Reference (7) 



ORIGINAL PAGE IS 
OF POOR QUALITY 



(b) Airfoil 


Figure 17. Concluded 


rteiGlNAt PAGE 



Figure 18. Comparison with Reference (7), calculated 
versus experimental . 

(a) Reference Figure = 2.0 
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QRieiNAL PAGE IS 
m POOR QUALITY, 


Experimental 

O Calculated 




;'oS5 QUAUTY 



1 1 1 p. 

0.2 0.4 0.6 0.8 1.0 

x/c 


Figure 19. Experimental chordwise pressure 
distributions from Reference (7). 
Reference Figure V-la, Vj =2.5 


i 
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— — Experimental 
O Calculated 



t 


I 


“1 — ~i r r“ 

0.4 0.6 0.8 1.0 

x/c 


Figure 21. Cotnparison with Reference (7), 
calculated versus experimental. 
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